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ABSTRACT 

An experimental aerodynamic evaluation of three 
compact radial turbine builds was performed. Two ro- 
tors which were 40-50 percent shorter in axial length 
than conventional state-of-the-art radial rotors were 
tested. A single nozzle design was used. One ro- 
tor was tested with the nozzle at two stagger angle 
settings. A second rotor was tested with the noz- 
zle in only the closed down setting. Experimental re- 
sults were compared to predicted results from a quasi- 
3D inviscid and boundary layer analysis, called MTSB 
(Meridl/Tsonic/Blayer). This analysis was used to pre- 
dict turbine performance. It has previously been cali- 
brated only for axial, not radial, turbomachinery. The 
predicted and measured efficiencies were compared at 
the design point for the three turbines. At the design 
points the analysis overpredicted the efficiency by less 
than 1.7 points. Comparisons were also made at off- 
design operating points. The results of these compar- 
isons showed the importance of an accurate clearance 
model for efficiency predictions and also that there are 
deficiencies in the incidence loss model used. 

INTRODUCTION 

Radial turbines have the advantage over axial tur- 
bines in their ability to extract high work per stage. 
They also may be attractive from a packaging stand- 
point. Unfortunately, the bulkiness and weight of ra- 
dial inflow turbines place them at a disadvantage in air- 
craft propulsion systems. To address this issue, a joint 
program was established between Pratt and Whitney 
Aircraft and NASA Lewis Research Center to signifi- 
cantly reduce the axial length, and thereby weight, of 
a radial-in-flow turbine rotor. 
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Under a Cooperative Endeavor Agreement, Pratt and 
Whitney designed two highly loaded compact radial- 
in-flow turbine rotors. Both parties analyzed the ro- 
tor flow fields with their respective design systems, 
and NASA Lewis conducted aerodynamic performance 
tests. 

The program goal was to reduce the rotor axial 
length by 40-50 percent and the weight by 20-30 per- 
cent while achieving equal or higher efficiency compared 
to conventional state-of-the-art radial rotors. The re- 
duction in axial length was accomplished primarily by 
shortening the exdueer portion of the blade, which nor- 
mally has the lowest aerodynamic loading. Carefully 
controlling the flow to avoid high losses in such a sig- 
nificantly shorter blade was made possible with the use 
of modern three- dimensional flow analysis codes. A 
nozzle and two rotors were designed to turbine require- 
ments suitable for a main propulsion engine in rotor- 
craft and/or regional commuter aircraft. The design 
approach and performance results for the compact rotor 
turbines are reported in reference 1. The experimental 
performance tests conducted in the Lewis Warm Tur- 
bine Facility 1 verified stage efficiencies equal or better 
than conventional rotors for this compact rotor concept. 

To further understand the aerodynamic charac- 
teristics of the compact rotor concept, a fully three- 
dimensional (3-D) viscous analysis and a quasi 3-D in- 
viscid/boundary layer analysis were made at the Lewis 
Research Center. The two analytical studies utilized 
entirely different computer codes and have different 
purposes. The viscous analysis code is a fully 3-D “av- 
erage passage” formulation that calculates the entire 
blade passage flow field and is used to analyze stator- 
rotor interaction effects of one or more stages. The 
results of that analysis for the compact radial turbine 
are reported in reference 3. The quasi-3D analysis sys- 
tem, which is the subject of this paper, is a rapid run- 
ning blade design and loss analysis tool. It calculates 



the inviscid flow field, the surface boundary layer and 
the blade losses, using empirical correlations, for a sin- 
gle blade row. The system of computer codes used for 
the inviscid analysis is referred to as MTSB and is de- 
scribed in reference 4. Additional details of MTSB and 
the analysis methodology are given in the Analysis Pro- 
cedure. 

The analysis of the compact turbine with MTSB 
was conducted to identify and quantify the loss sources, 
determine the differences between the three configura- 
tions and indicate possible areas of performance im- 
provement. In addition, comparison of the calculated 
stage performance with measurements were used to 
evaluate the MTSB loss system for radial turbines. The 
loss correlations in MTSB are heavily based on axial 
turbine data and predict the performance of those tur- 
bines very well, reference 4. The application of this 
procedure to the design of a highly contoured radial 
turbine is described in reference 5. However, test data 
were not available to assess the prediction’s accuracy. 
The comparison of the MTSB performance prediction 
results with experimental data for a conventional radial 
turbine is discussed in reference 6. 

This report describes the application and results of 
the quasi- three dimensional analysis of the three com- 
pact radial turbines tested. The calculated loss con- 
stituents are tabulated and overall performance com- 
pared to the experimental results. Comparison of the 
analysis with experimental data is made at the design 
speed, weight flow, and pressure ratio. In addition, 
MTSB was used to calculate the off-design losses at 
design speed for one turbine configuration and the re- 
sults compared to test values. Stator and rotor blade 
surface static pressures were also calculated and com- 
pared for the inviscid and viscous analyses. Finally, 
improvements to the MTSB loss models applicable to 
radial turbines are discussed. 

TURBINE DESCRIPTION 

Figure la shows the two compact radial turbine 
rotors along with a conventional rotor. The compact 
turbines are approximately 42 and 52 percent shorter in 
axial length compared to the conventional rotor. Each 
rotor had 14 blades. The backface of the rotor is scal- 
loped, i.e., the backface disk does not extend all teh way 
to the leading edge of the blade for weight and stress 
considerations. With a scalloped hub flow leakage can 
occur in the backface region as well as in the blade 
tip shroud region. The stator, which had 36 vanes, is 
shown in figure lb. The flow enters radially through the 
inlet vanes and is turned tangentially approximately 67 


degrees with the design stator. 

Three turbine stage configurations were evaluated 
in the experiment and analyzed with MTSB. They are; 
(1) Rotor I with the design inlet vanes; (2) Rotor I with 
the stator vanes closed down 1.125 degrees; (3), Rotor 
II with the closed down vanes. The three configurations 
are referred to in the text as Rotor I Stator I, Rotor I 
Stator II, and Rotor II Stator II, respectively. Rotor II 
was designed to be physically smaller than rotor I and 
therefore had a lower blade speed for the same rotative 
speed. However, both rotors were scaled up to a tip 
diameter of .367 meters (14.59 inches) for the experi- 
mental evaluation and the rpm of rotor II decreased to 
match the design equivalent conditions. The test con- 
ditions for the three configurations are given in table I. 
Additional design information is contained in reference 
1 . 

ANALYSIS PROCEDURE 

The analysis used the quasi-3D procedure devel- 
oped for axial turbines by Boyle, Haas, and Katsanis^ 4 ). 
The quasi-3D flow analysis was done by first obtain- 
ing a hub-toshroud midchannel flow solution using a 
modified version for the MERIDL code developed by 
Katsanis and McNally^ 7 ). The results were then used 
to obtain a series of blade-to-blade solutions using the 
TSONIC code developed by Katsanis^8). The flow so- 
lutions were then used as input to a modified version 
for the BLAYER boundary layer analysis of McNally^ 9 ). 
The results of the boundary layer analyses were used 
along with correlations for incidence, tip clearance and 
secondary flow loss to obtain the predicted efficiency of 
the stage. 

Civinskas and Povinelli^ 5 ) made additions to the 
procedure to extend the analysis to radial turbines. 
They added Wiesner’s^ 10 ) optimum incidence model, 
and used the radial turbine clearance model of Futral 
and Holeski^ 11 ) to predict efficiencies in a radial turbine 
designed to avoid separation. 

The clearance model used by Civinskas and 
Povinelli did not account for any additional loss due 
to a backface scallop. A radial clearance loss model, 
similar to the one proposed by Rodgers^ 12 ), was used 
herein in the efficiency predictions. In this clearance 
loss model the change in stage efficiency, Atjcl, as a 
fraction of the zero clearance efficiency, rj 0 , is given as 
a function of three rotor clearances. This model is: 

= L a Ca + LrCr + L b C b (1) 

* 7 ° 

C a is the axial clearance, C r is the radial clearance, 
and Cb is the backface clearance. Both C a and Cb are 
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found by dividing the actual axial and backface clear- 
ances by the inlet span. C r is the ratio of the radial 
clearance to the span at the exit. The coefficients L a 
and L r were given as 0.15 and 1.6 respectively by Fu- 
tral and Holeski^ 11 ) for a radial turbine with no scal- 
lop. Rogers^ 12 ) stated that the equivalent shroud axial 
clearance equals the front shroud plus one half the back 
shroud clearance. This approach, if used in the present 
designs, would result in a very small efficiency penalty 
due to the presence of the scallop. The efficiency loss 
due to the scallop would be less than 15% of the to- 
tal efficiency loss due to clearance. The coefficient L\, 
was instead taken as 0.8, resulting in the efficiency loss 
due to the scallop clearance being nearly 60% of the 
total clearance loss. In axial turbines the clearance 
loss is often taken as a function of rotor tip reaction 
(Glassman^ 13 )). No information was found in the open 
literature regarding the effect of tip loading on rotor 
clearance loss for radial turbines, and so none was used 
in the present analysis. 

Secondary loss correlations such as Morris and 
Hoare^ 14 ), and Dunham^ 15 ) account for both the loss 
due to an inlet boundary layer and the loss due to end- 
wall boundary layer. The analysis used herein explicitly 
accounts for the endwall boundary layer loss, and uses 
a correlation only for the loss due to the inlet bound- 
ary layer. The loss correlation used is the same as that 
used in reference 4, and is: 


y S EC = 0.44Z 


COSOL 2 / Cw^n 

cos (3 1 '/r ' c ' 


( 2 ) 


As observed by Civinskas and Povinelli^ 5 ) one for 
the deficiencies of the analysis procedure of reference 4 
is the treatment of separated flow. In this procedure, 
when the predicted velocity distribution resulted in sep- 
aration, the velocity distribution was smoothed suffi- 
ciently so that no separation occurred. Losses were cal- 
culated using the boundary layer parameters calculated 
using the smoothed velocities. This resulted in higher 
predicted losses for cases without separation than for 
cases with massive separation. To remedy this situa- 
tic^jj the boundary layer analysis procedure was modi- 
fied. For the analysis used herein, whenever separation 
occurred, the boundary layer parameters were held con- 
stant from a point just upstream of separation to a lo- 
cation where a favorable pressure gradient indicated no 
separation. The parameters which were held constant 
in the separated region were the boundary layer thick- 
nesses and form factors. 

The MTSB performance prediction analysis is per- 
formed using experimentally measured values of tem- 


perature, pressure, rotative speed, and massflow as in- 
put. To obtain the MTSB stage solution the blade rows 
are run sequentially. The stator performance and exit 
swirl are first determined. For transonic flows a check 
is made to insure that the ratio of the stator exit static 
pressure to the inlet total pressure is consistent with the 
experimental value. The stator exit conditions are then 
used as input to the rotor analysis. If the calculated 
stage total pressure ratio differs from the experimental 
one, conditions are modified until the desired pressure 
ratio is achieved. The primary convergence criteria for 
each blade row is that the trailing edge Kutta condition 
is satisfied along the entire span. Losses affect the flow 
solution for each blade row, and are updated during the 
iteration for the blade row exit swirl. 

RESULTS 

Figure 2 shows the computational mesh used for 
analyzing Rotor I Stator I hub-to-shroud flow. Similar 
meshes were used to analyze the other two stages. 

SURFACE VELOCITIES 

Blade surface velocities are shown in figures 3 
through 5 depicting Rotor I Stator I, Rotor I Stator II, 
and Rotor II Stator II configurations respectively. The 
velocities shown are at the hub, mean and tip regions 
of each blade. 

Stator I accelerates the flow to a maximum value 
of velocity V/Vcr=1.2 at approximately 70 percent of 
the blade meridional distance. The maximum veloci- 
ties occur at the hub and tip regions. Stator II has a 
maximum velocity ratio of 1.3 due to the closed down 
nature of the flow. For both rotors the stator velocities 
are similar with Rotor II resulting in a slightly higher 
peak stator velocity ratio. 

The rotor surface velocities show that after the flow 
accelerates at the leading edge of the blade, diffusion 
takes place. The inlet velocities are higher for Rotor 
II than for Rotor I. However, the highest surface ve- 
locity occurs on Rotor I at the tip of the blade in the 
exducer region. Analysis showed small areas of separa- 
tion in the rotors for all three configurations. Rotor I 
Stator I showed a separation region starting at approx- 
imately 37 % of the meridional distance. Due to the 
favorable pressure gradients farther downstream, the 
analysis predicted reattachment at approximately 42% 
of the meridional distance. The Rotor I Stator II con- 
figuration showed a similar separation region near the 
hub. For this configuration the flow also separated in 
the tip region near the trailing edge. Rotor II had hub 
region separation between 50% to 55% of the merid- 
ional distance, and no separation in the tip region. 
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BLADE STATIC PRESSURES 

Figures 6 through 8 show the blade loading com- 
puted for the three configurations. Surface static pres- 
sure contours are shown normalized with inlet total 
pressure, (Pf o). The stator blades shown have min- 
imum static-to-total pressure ratios between .35 and 
.45, The minimum static pressure occurs on the suc- 
tion surface at approximately 70 percent of the chord. 
High adverse pressure gradients occur at the exit of the 
stator on the suction side of the vanes. 

The static pressure distributions are determined 
from the inviscid velocity distributions. Because the 
total pressure is nearly constant, the pressures stream- 
wise are nearly the same as those calculated from the 
isentropic relations. 

In contrast with stators, the rotor has a varying 
relative total pressure change due to the change in ra- 
dius with meridional distance. The rotor pressure con- 
tours on the suction side show that the flow at the inlet 
has high pressure gradients. These gradients are due to 
the change in relative total pressure, as well as the de- 
creasing static pressure with the radius. Comparing 
the suction and pressure surface pressure distributions 
shows that the rotor blades have highest loading in the 
shroud region. The suction side of Rotor II contours 
show low pressure areas near the inducer. These areas 
are located near the back face of the rotor where the 
blade is scalloped. The low pressure shown in the Ro- 
tor II inducer backface area is unique to Rotor II. All 
rotor configurations show a low pressure area at the tip 
region in the bend of the blade. These regions of low 
pressure on the suction side of the blade can result in 
large leakage of flow from the pressure side. 

A comparison was made between the calculated 
MTSB surface static pressures made and those calcu- 
lated by Heidmann^ 3 ) using a viscous fully 3-D “average 
passage” code. The stator comparisons are shown in 
figure 9, and the rotor comparisons are shown in figure 
10. The comparison was made for the Rotor I Stator 
I configuration. These figure shows three streamlines 
near the hub, at midspan, and near the tip. Figure 9 
shows that the fully 3D viscous solution predicts the 
stator minimum suction side pressure to occur further 
downstream than the MTSB analysis. However, both 
analyses show nearly the same minimum pressure level 
at each streamline. Figure 10 shows similar behavior for 
the rotor. The differences may be due to differences in 
the loss distribution. An empirical distribution of losses 
was used in the MTSB analysis. The fully 3D viscous 
analysis accounted for all clearances. While the MTSB 
analysis predicted separation for this rotor configura- 


tion, the fully 3D solution did not. A slight variation 
in exit pressures used to run the two cases caused the 
difference in the loading at the trailing edges of both 
stator and rotor. Overall the agreement with the two 
calculations is reasonable. 

LOSSES 

Figure 11 compares the predicted and experimen- 
tal loss for the three configurations at design conditions. 
Also shown is a breakdown of the predicted loss. The 
first two configurations have nearly the same test effi- 
ciency. Configuration three had about a one-half point 
lower efficiency. The analysis predicted a variation of 
less 0.6 points in efficiency among the three configura- 
tions. The analysis underpredicts the experimental loss 
by less than 1.7 points. This was similar to the results 
of reference 4 where the analysis predicted the design 
point efficiency within 1.2 points for nine axial turbine 
stages. 

The loss breakdown shows that clearance loss is 
the major predicted loss mechanism. As discussed pre- 
viously, the clearance loss model adopted was expected 
to be conservative with respect to the loss prediction. 
It was felt that the model was more likely to overpredict 
the loss due to clearance, than to underpredict this loss. 
However, the data show higher losses than predicted. 
This indicates that if the clearance loss model is con- 
servative, it is not unduly conservative. 

The stator loss is the next largest loss component 
at the design. This loss is the sum of profile, endwall, 
and secondary flow losses. Approximately 43% of the 
stator loss was profile loss. The endwall loss, which 
accounted for boundary layer growth up to the rotor 
leading edge was 50% of the total loss. The remainder 
is due to secondary loss, which is small because of the 
large span at the stator inlet. 

The rotor profile loss is also a significant loss. The 
profile loss is calculated from the trailing edge blockage 
and boundary layer growth. It is not expected to be 
a major source of uncertainty in the loss calculation. 
The endwall and secondary loss together are smaller 
than the predicted rotor profile loss. The rotor end- 
wall loss is small because the exit span is large relative 
to the hub endwall momentum thickness. The rotor 
secondary loss is small because of the acceleration in 
the vaneless space. While the endwall and secondary 
losses have higher uncertainties than the profile loss, 
their uncertainties are not sufficient to account for the 
differences between the measures and predicted stage 
losses. 

It is somewhat surprising to see a significant inci- 
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dence loss at the design point. The reason for relatively 
high loss at design is that the incidence loss model pre- 
dicts that the optimum incidence angle occurs when 
the ratio of absolute tangential velocity to wheel speed, 
( V u /U ), is less than 1. Other incidence loss models for 
radial turbines also have this characteristic. However, 
at the design point V u /U was predicted to be greater 
than 1. This was true for all three configurations. The 
measured values of the stator pressure ratio confirmed 
this. Figure 12 shows the predicted and measured effi- 
ciency as a function of overall pressure ratio at the de- 
sign speed for the Rotor II Stator II configuration. The 
comparisons were similar for all three configurations. 
As the stage pressure ratio decreased, the stator pres- 
sure ratio decreased. Consequently, V u /U decreased, 
and the predicted incidence loss decreased. This re- 
sulted in increased predicted efficiency at lower pressure 
ratios. The experimental efficiency remained fairly con- 
stant over the range of pressure ratios. These results 
indicate a deficiency in the incidence loss model. A 
partial explanation for the inaccuracy of the incidence 
model might be the use of elliptical leading edge for 
both the stator and rotor blades. All of the incidence 
loss models are based on circular leading edges. 

SUMMARY AND CONCLUSIONS 

Three compact radial turbines were analyzed using 
MTSB, a quasi- 3D inviscid and boundary layer anal- 
ysis. The results of the analysis were compared with 
experimental data. At design conditions MTSB pre- 
dicted the test efficiencies within 1.7 points for all three 
configurations. The calculated loss sources suggest that 
designs, such as these tested, which have optimum effi- 
ciencies at high values of V u /U , can result in efficiency 
gains of up to two points due to minimization of the 
incidence loss. 

Comparisons of blade surface static pressures cal- 
culated using MTSB agreed reasonably well with those 
calculated using a fully three-dimensional viscous code. 
However, MTSB predicted some small regions of sep- 
aration in the rotors, which were not predicted using 
the viscous analysis. The turbine test program did not 
provide information to confirm either flow situation. 

The use of MTSB to predict the off design perfor- 
mance of these compact radial turbines was not entirely 
successful. As the stage pressure ratio was decreased 
along a constant design speed line, MTSB predicted an 
increase in efficiency. However, the experimental data 
show nearly constant efficiencies over the same range 
of pressure ratios. The cause of this discrepency was 
the incidence model being used. This discrepency indi- 


cates the need for an improved incidence model. The 
clearance loss model used for the radial turbine perfor- 
mance calculation gave reasonable results. The impor- 
tance of the clearance loss to the overall loss shows the 
importance of an accurate clearance loss model. A pre- 
cise model would include more factors than were used 
in the simplified clearance model used in this analysis. 
Overall the changes made to MTSB to extend its use to 
radial turbines has enhanced its usefulness in analyzing 
the advanced compact radial turbines. 
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Rotor 1 Stator 1 

Rotor 1 Stator II 

Rotor II Stator II 

ro,R 

860.000 

860.000 

860.000 

P'0, psia 

48.400 

48.400 

45.500 

Spec. Work, Btu/lbm. 

69.400 

66.030 

67.870 

Work Factor 

1.091 

1.048 

1.195 

D tip, in. 

14.459 

14.459 

14.459 

Rotor Tip Speed, f/s 

1,263.200 

1.257.300 

1.194.000 

Rotative Speed, RPM 

20,022.500 

19.928.500 

18,926.400 

Mass Flow Rate, Ib/s 

6.295 

5.940 

5.630 

Power, HP 

617.970 

554.780 

540.510 

P'in/P’out 

5.000 

5.000 

5.000 


Table I.- Turbine test conditions 


6 
































Figure la.- Compact and conventional size turbine 

rotors 


Figure lb.- Stator vanes 
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Figure 2.- MTSB computational mesh 
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Stator Surface Velocities 


Rotor Surface Velocities 
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Figure 3.- Rotor I Stator I calculated blade surface velocities 


Stator Surface Velocities 


Rotor Surface Velocities 




Figure 4.- Rotor I Stator II calculated blade surface velocities 
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Figure 5.- Rotor II Stator II calculated blade surface velocities 
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Figure 6.- Calculated blade loading, surface static pressure contours for Rotor I Stator I 
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Figure 7.- Calculated blade loading, surface static pressure 


contours for Rotor I Stator II 
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Figure 8.- Calculated blade loading, 
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static pressure contours for Rotor II Stator II 
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Figure 9.- Stator surface static pressures for MTSB 

and viscous 3D predictions, Rotor I Stator I. 
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Figure 10.- Rotor surface static pressures for MTSB 

and viscous 3D predictions, Rotor I Stator I 
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Figure 11.- Comparison of experimental and computed losses for the three turbine stages 



Figure 12.- Predicted and measured efficiencies for Rotor II Stator II 
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